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ABSTRACT 

We present in this paper a detailed analysis of the effect of environment on the star-formation activity of 
' galaxies within the Early Data Release (EDR) of the Sloan Digital Sky Survey (SDSS). We have used the Ha 

, emission line to derive the star-formation rate (SFR) for each galaxy within a volume-limited sample of 8598 

galaxies with 0.05 < z < 0.095 and M(r*) < -20.45. We find that the SFR of galaxies is strongly correlated with 
the local (projected) galaxy density and thus we present here the density-SFR relation that is analogous to the 
, density-morphology relation. The effect of density on the SFR of galaxies is seen in three ways. First, the overall 

distribution of SFRs is shifted to lower values in dense environments compared with the field population. Second, 
the effect is most noticeable for the strongly star-forming galaxies (Ha EW > 5A) in the 75 th percentile of the 
SFR distribution. Third, there is a "break" (or characteristic density) in the density-SFR relation at a local galaxy 
density of ~ 1 /i^Mpc" 2 . To understand this break further, we have studied the SFR of galaxies as a function of 
—i , clustercentric radius from 17 clusters and groups objectively selected from the SDSS EDR data. The distribution 

■ of SFRs of cluster galaxies begins to change, compared with the field population, at a clustercentric radius of 3-4 

virial radii (at the > lcr statistical significance), which is consistent with the characteristic break in density that we 
0^ ■ observe in the density-SFR relation. This effect with clustercentric radius is again most noticeable for the most 

strongly star-forming galaxies. 

Our tests suggest that the density-morphology relation alone is unlikely to explain the density-SFR relation 
we observe. For example, we have used the (inverse) concentration index of SDSS galaxies to classify late-type 
galaxies and show that the distribution of the star-forming (EW Ha > 5 A) late-type galaxies is different in dense 
regions (within 2 virial radii) compared with similar galaxies in the field. However, at present, we are unable to 
make definitive statements about the independence of the density-morphology and density-SFR relation. 

We have tested our work against potential systematic uncertainties including stellar absorption, reddening, 
SDSS survey strategy, SDSS analysis pipelines and aperture bias. Our observations are in qualitative agreement 
with recent simulations of hierarchical galaxy formation that predict a decrease in the SFR of galaxies within the 
virial radius. Our results are in agreement with recent 2dF Galaxy Redshift Survey results as well as consistent 
with previous observations of a decrease in the SFR of galaxies in the cores of distant clusters. Taken all together, 
these works demonstrate that the decrease in SFR of galaxies in dense environments is a universal phenomenon 
•rH . over a wide range in density (from 0.08 to lO/j^lMpc" 2 ) and redshift (out to z — 0.5). 

Subject headings: galaxies: clusters: general - galaxies: evolution - stars: formation - galaxies: stellar content - 
C$ ' surveys 
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1. INTRODUCTION 

In this paper, we investigate the relation between environ- 
ment and star formation rate (SFR) of nearby galaxies. Spec- 
troscopic studies of distant clusters of galaxies (z > 0.2) have 
already established a clear connection between these two phys- 
ical properties, in that the SFR of galaxies in the cores of dis- 
tant clusters are significantly lower than those observed in the 
field at the same redshift (Balogh et al. 1997; Hashimoto et al. 
1998; Poggianti et al. 1999; Couch et al. 2001; Postman et al. 
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2001). In particular, Balogh et al. (1998) found that the SFRs of 
cluster galaxies were lower relative to field galaxies of similar 
bulge-to-disk ratio, physical disk size, and luminosity, which 
suggests that the observed decrease in the star-formation may 
not be fully explained by the density-morphology (Dressier 
1980) or radius-morphology (Whitmore et al. 1993) relations 
(see also Poggianti et al. (1999); Couch et al. (2001); Balogh 
et al. (2002a)). This conclusion has also been suggested by 
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Hashimoto et al. (1998), who used the Las Campanas Redshift 
Survey (LCRS) to study the relationship between local galaxy 
density, SFR, and galaxy concentration index, and found that 
cluster galaxies have a reduced level of SFR compared with the 
field, regardless of their concentration index. All these stud- 
ies provide circumstantial evidence that galaxies are undergo- 
ing a physical transformation as they enter dense environments, 
and that the timescale for this transformation is different for the 
morphological and star formation properties. 

In models of hierarchical galaxy formation (Kauffmann et al. 
1993; Somerville & Primack 1999; Cole et al. 2000), it is as- 
sumed that galaxies accreted into larger halos (e.g., they fall 
into groups or clusters) have their hot gas reservoir removed, 
which results in a gradual decline of the star-formation activity 
of galaxies in dense environments. This simple prescription is 
able to successfully match some of the general observed trends 
(Baugh et al. 1996; Diaferio et al. 2001; Balogh et al. 2000). 
At present, however, these models do not explore the physical 
details of how the galaxies lose their gas (see Bekki, Shioya & 
Couch 2002). Many possible mechanisms are proposed in the 
literature, such as ram pressure stripping of the cold gas due to 
interaction with the intracluster medium (Gunn & Gott 1972; 
Quilis et al. 2000), galaxy harassment from high velocity en- 
counters with other galaxies (Moore et al. 1999), tidal disrup- 
tion (Byrd & Valtonen 1990) or galaxy evolution via mergers 
and close encounters in in-falling poor groups (Zabludoff et al. 
1996; Zabludoff & Mulchaey 1998). Proper treatment of these 
effects will likely affect the details of the predicted relationships 
between galaxy properties and environment, but the relative im- 
portance of the proposed evolutionary mechanisms can only be 
determined with better data than has been available so far. 

As a first step in addressing these issues, we use the Early 
Data Release (EDR; Stoughton et al. 2002) of the Sloan Digital 
Sky Survey (SDSS; see York et al. 2000) to further understand 
the relation between environment and the star-formation rate 
of galaxies. The EDR is an excellent database for such a study 
because it provides a large, homogeneous sample of galaxies, 
with medium resolution, flux-calibrated spectra and 5 passband 
photometry, over a wide range of environments. Such data al- 
low for the computation of the Ha and [Oil] equivalent widths 
(EW) and SFR for all galaxies in a self-consistent way. Our 
work expands upon earlier cluster studies, as discussed above, 
in four key ways. First, we probe the rarefied field and poor 
group regimes, which include the most common galactic envi- 
ronments in the Universe. Second, we focus on the low redshift 
universe (z < 0. 1), which provides a zero-point for the high red- 
shift studies. Third, we are able to follow the star-formation of 
galaxies continuously from the cores of rich clusters and poor 
groups into the rarefied field. Finally, we can accurately quan- 
tify the local galaxy density in a uniform manner that is not 
subject to statistical background or luminosity corrections. Our 
work is complementary to earlier field studies using the 2dF 
Galaxy Redshift Survey (2dFGRS, see Lewis et al. 2002) and 
the LCRS (Hashimoto et al. 1998) and explores the density- 
SFR relation in greater detail than these previous works. 

In Section 2, we present a description of the SDSS EDR data, 
the sample selection, and our SFR measurements. Our results 
are presented in Section 3. In Section 4, we discuss possible 
systematic biases, the density-morphology relation, and mod- 
els of hierarchical structure formation. We also compare our 
work with previous studies, e.g., the 2dFGRS. In Section 5, we 
present our conclusions. In the Appendix, we discuss the de- 



tails of selecting clusters and groups from the SDSS data. We 
also consider the effects of systematic biases that may influence 
our results, including the internal reddening of galaxies, sample 
selection, data analysis pipelines, and luminosity and aperture 
biases. Throughout this paper, we use H = 75kms _1 Mpc _1 , 
f2 m = 0.3 and 17a = 0.7 unless otherwise stated. 

2. DATA 

2.1. Early Data Release of the Sloan Digital Sky Survey 

The Sloan Digital Sky Survey (http://www.sdss.org) is a 5 
passband (u*, g*, r*, ;*, z*) imaging and medium-resolution 
(R ~ 1800) spectroscopic survey of the Northern Galactic 
Hemisphere (see York et al. 2000, for details). For technical de- 
tails regarding the SDSS imaging survey, the reader is referred 
to Gunn et al. (1998), Smith et al. (2002), Pier et al. (2002) and 
Hogg et al. (2002). In June 2001, the SDSS publicly released 
photometric and spectroscopic data comprise of nearly 50,000 
spectra of galaxies, stars and QSOs over 460 deg 2 of sky. This 
release is known as the Early Data Release (EDR) and is fully 
described in Stoughton et al. (2002). 

One of the unique features of the SDSS is the set of sophis- 
ticated data analysis software pipelines (see Lupton et al. 2001, 
Stoughton et al. 2002 for details) used to reduce the raw im- 
ages and spectra into large catalogs of sources. Briefly, for 
each flux-calibrated SDSS spectrum, a redshift is determined 
both from the absorption lines (via cross-correlation; Heavens 
1993) and the emission lines (via a wavelet-based peak-finding 
algorithm; Frieman et al. in prep). Once the redshift is known, 
the spectroscopic pipeline estimates the continuum emission 
at each pixel using the median value from a sliding box of 
100 pixels (~ 100A) centered on that pixel. Emission and ab- 
sorption lines are measured through the fitting of a Gaussian, 
above the best-fit continuum, at the redshifted rest-wavelength 
of the lines. In order to accommodate line blending, the SDSS 
pipeline fits multiple Gaussians in the case of the Ha and [Nil] 
doublet and the [OIII] doublet. Thus, for all the major emis- 
sion/absorption lines in the galaxy spectra, the spectroscopic 
pipeline provides an estimate of the equivalent width (EW), the 
continuum level (at center of the line), a line identification (e.g., 
Ha ), a goodness-of-fit (% 2 ), and the height and sigma of the 
fitted Gaussian (and the associated statistical errors on all these 
quantities). These quantities are measured for all the major stel- 
lar emission and absorption lines regardless of their formal de- 
tection. In the Appendix, we present preliminary quality assur- 
ance tests of these emission line measurements which show that 
the Ha and [Oil] fluxes and EWs are robust. 

When converting magnitude to flux, we have ignored the 
asinh magnitudes of Lupton et al. (1999) and treated the 
SDSS-EDR Petrosian magnitudes as traditional AB magni- 
tudes (Fukugita et al. 1996). At the bright magnitudes used 
herein, the difference between these two magnitude systems is 
less than 1%, i.e., less than the photometric accuracy of the 
SDSS (see Stoughton et al. 2002). 

2.2. Sample Selection 

The reader is referred to Strauss et al. (2002) for a detailed 
description of the spectroscopic target selection for the SDSS 
main galaxy survey. For the analyses presented in this paper, we 
have used the EDR data discussed above but with the follow- 
ing additional selection criteria. We begin with all objects that 
have been spectroscopically confirmed as galaxies and have a 
redshift confidence of > 0.7 (see Heavens 1993). These criteria 
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result in a sample of 41,622 galaxies. We then reject 37 ob- 
jects that have certain redshift warning flags set in the database, 
i.e., we exclude spectra with the Z_WARNING_NO_BLUE (no 
blue side of the spectrum) and Z_WARNING_NO_RED (no 
red side of the spectrum) flags (see Stoughton et al. 2002). We 
also reject 1231 galaxies because they are duplicates. We have 
excluded another 37 galaxies that have z* > 22.83, as these 
are likely spurious detections in the SDSS (see Table 21 of 
Stoughton et al., 2002), and, at this magnitude, the effect of 
the asinh magnitudes starts to become noticeable. This selec- 
tion leaves us with 40,317 galaxies. Finally, we exclude 4014 
galaxies for which the Ha emission line could not have been 
measured and 3452 galaxies for which the [Oil] emission could 
not have been measured. This is because these emission lines 
fall in masked regions of the spectrum, due to cosmic ray hits, 
cosmetic defects in the CCDs or missing data. 

We define a volume-limited sample by further restricting the 
sample to galaxies in the redshift range 0.05 < z < 0.095 
and more luminous than M{r*) = -20.45 (k-corrected, for 
H Q = 75 km s -1 Mpc -1 ). This magnitude limit corresponds to 
~ M*(r*) + 1, assuming M*» = -20.8 + 51og/i from Blanton 
et al. (2001). The lower redshift limit is imposed to minimize 
aperture bias (see the Appendix), while the upper limit is the 
redshift where our luminosity limit equals the magnitude limit 
of the SDSS (r* = 17.7; Strauss et al. 2002). These cuts leave 
us with a volume-limited sample 8598 galaxies for our anal- 
yses and includes all types of galaxies, not just emission line 
objects. 

2.3. Measuring the Star-Formation Rate in Galaxies 

The star formation rates of galaxies can be indirectly deter- 
mined in various ways, as reviewed by Kennicutt (1998a). The 
most successful methods rely on correlations of the star for- 
mation rate with measurements of the far-infrared luminosity 
(which arises from dust heated by star formation), the radio lu- 
minosity (which results from synchroton emission associated 
with supernovae) and indicators that are sensitive to the ion- 
izing flux from massive stars. The last category includes mea- 
surements of the ultraviolet continuum and the fluxes of nebular 
emission lines. Each of these techniques is subject to different 
biases and calibration uncertainties, but does give consistent es- 
timates for "normal" galaxies when these effects are accounted 
for (Chariot et al. 2001; Hopkins et al. 2001; Sullivan et al. 
2001). In this paper we focus on SFRs estimated from the Ha 
line. We also present the [Oil] EW measurements to facilitate 
comparison with high redshift studies of galaxies. Both lines 
are sensitive to the metallicity and ionization levels of the gas, 
though the problem is more significant for [Oil]. Dust extinc- 
tion is the largest source of uncertainty, as both the effective 
optical depth and dust geometry must be considered (Chariot & 
Longhetti 2001). 

In Figure 1, we show the distribution of rest-frame Ha and 
[Oil] EWs, for the sample of galaxies discussed above. There 
is a small tail of negative EWs caused by statistical errors and, 
in the case of Ha , stellar absorption. Typically, stellar absorp- 
tion is expected to account for < 5 A of the Ha EW (Kennicutt 
1992; Chariot et al. 2001) and is therefore a small contribution 
to the majority of the strong emission line galaxy spectra (see 
the Appendix for more details). The distributions of [Oil] and 
Ha EWs are highly asymmetrical, skewed to low values but 
with a long tail to large, positive values. The observed mean of 
the SDSS [Oil] distribution (see Figure 1) is half the value of 



the observed mean for the field galaxy sample of the CNOC1 
distant cluster studies (Balogh et al. 1997). This difference may 
be due to the differences in the apertures used by the SDSS and 
CNOC (although at z=0. 1 the SDSS aperture corresponds to 5.5 
kpc, while at z=0.3 the CNOC aperture corresponds to 6.7 kpc), 
or the evolution in the SFR of field galaxies (Lilly et al. 1996). 

We have derived SFRs for these galaxies using the Ha flux 
(as computed from the pipeline output) and the theoretical rela- 
tion from Kennicutt (1998b), 

SFR(M Q /yv) = 7.9 x 10 -42 L(Ha), (1) 

where L(Ha) is the observed luminosity in the Ha line (in erg 
s" 1 ). This relation is valid for so-called Case B (optically thin) 
recombination, with a Salpeter (1955) initial mass function. We 
make no correction for negative SFRs that are caused by the 
small tail of negative Ha EWs in Figure 1 . Although nonphysi- 
cal, the negative SFR of these galaxies are included when com- 
puting the median, 25'*, and 75 th percentiles of the SFR distri- 
butions because removing them, or setting them to zero, would 
artificially skew the distributions. Alternatively, we could add 
a constant amount of stellar absorption to our sample but this 
would still leave a small negative tail, due to statistical errors, 
and the exact amount of stellar absorption to add remains un- 
clear (see Appendix). We note here that the SFRs presented 
in this paper are not corrected for the 3 arcsecond SDSS fiber 
aperture and are therefore, systematically lower, by a factor of 
~ 5, compared to total SFRs derived from the radio or by inte- 
grating the light from the whole galaxy (see Hopkins et al., in 
prep) 

We then correct the SFR for SFR-dependent reddening us- 
ing the empirical formulae of Hopkins et al. (2001). We have 
also computed the normalized SFR (SFRN; see Postman et al. 
2001), which is the observed SFR divided by the z* luminos- 
ity of each galaxy in our sample, i.e., SFRN = SFR/L^, B where 
L^, B is the k-corrected (< 0.2 magnitudes) AB z*-band lumi- 
nosity (in erg s" 1 A" 1 at the effective wavelength of the z* fil- 
ter, 9049 A) computed via the prescription outlined in Fukugita 
et al. (1996). To make the units more intuitive, we have re- 
normalized the z* luminosities by the characteristic z* AB 
luminosity from Blanton et al. (2001), i.e., L*(z*) = 6.34 x 
10 41 ergs _1 A" 1 . As discussed in Postman et al. (2001), such 
near-IR luminosities are less affected by recent star-formation 
and k-corrections and are more representative of the underly- 
ing stellar mass of the galaxy (at a given redshift) than are lu- 
minosities measured in bluer bands. The SFRN is therefore 
approximately the SFR per unit stellar luminosity, which helps 
to remove any large luminosity (or mass) bias from our results. 
We have also investigated the z-band luminosity function (LF) 
as a function of galaxy density and find, for the narrow range 
of bright z-band luminosities covered in this study, that the LF 
remains unchanged (within the errors) as a function of local 
galaxy density. This demonstrates that our results are unaf- 
fected by possible changes in the luminosity function of galax- 
ies with local galaxy density. 

3. RESULTS 

In this section, we present the results of two analyses of 
the relation between environment and the star-formation rate 
of galaxies. The first is an investigation of how the distribu- 
tion of galaxy SFRs changes as a function of local (projected) 
density. This analysis spans over two orders of magnitude in 
density, from the rarefied field (0.08 Mpc" 2 ) to the densest 
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environments (lO/z^Mpc -2 ). The second analysis focuses only 
on virialized systems with velocity dispersions ranging from 
~ 200 kms" 1 to 1000 kms" 1 and investigates the SFR of mem- 
ber galaxies as a function of clustercentric radius. 

3.1. SFR as a Function of Local Galaxy Density 

We investigate the SFR of galaxies as a function of their 
local galaxy density. We can parameterize the local environ- 
ment of each galaxy using the projected density of galaxies in 
a manner similar to that employed by Dressier (1980). This 
involves calculating the projected distance to the 10'* nearest 
spectroscopically-observed neighbor, and then converting that 
into a surface density (hj 2 Mpc~ 2 ). As discussed in Section 2.2, 
our limit is M(r*) = -20.45 (h = 0.75), which is close to the 
limit used by Dressier (M(V) = -20.4, h = 0.5). Our method- 
ology does differ from that used by Dressier (1980) in that 
we use a circular search aperture (instead of a rectangle) and 
spectroscopically eliminate foreground and background galax- 
ies so no further background correction is necessary. Thus 
our tenth nearest neighbor will be farther away than obtained 
via Dressler's definition. Therefore, our effective "smoothing" 
scale is larger, and we may not be as sensitive as Dressier (1980) 
to the densest environments, i.e., cluster cores. 

When calculating the local density, edge effects are particu- 
larly important because the EDR geometry is long strips of data 
that, at the median redshift of our sample, subtend only 14.8 
Mpc on the sky. To account for this, we remove any galaxy 
in which less than 75% of the circular area including the 10'* 
nearest neighbor lies within the EDR boundary. This correc- 
tions removes 1972 galaxies, leaving us with 6626 galaxies for 
our analysis. We do not make an area correction to our density 
estimates because we do not know if the 10'* nearest neighbor 
we have observed is, in fact, the true 10'* nearest neighbor. If 
it is, then correcting the density estimates for the area outside 
the survey boundary would bias the density estimates of those 
galaxies. Instead, we accept the fact that ~ 25% of galaxies 
near the edge of the survey have an incorrect local density mea- 
surement but that they are not systematically biased. We have 
checked the effect of changing our tolerance to accepting galax- 
ies near the edge of the survey, from rejecting all galaxies where 
the 10'* nearest neighbor is closer than any survey boundary, to 
including all galaxies regardless of edge effects, and find our 
results are robust. 

In Figure 2, we show the distribution of local galaxy den- 
sities for the 6626 galaxies that satisfy this constraint. These 
local galaxy density measurements are significantly lower than 
those presented by Dressier (1980) and Dressier et al. (1997). 
These differences are partly due to the cosmologies used, the 
definitions of the field and search area (see above), differences 
in the edge corrections, and probable errors in the magnitudes 
used by Dressier (1980), which have uncertainties of at least 0.5 
magnitudes. However, the most dominant effect is that we are 
probing lower density regions than Dressier (1980). 

In Figure 3, we present the Ha and [Oil] EW distributions, 
as quantified by the median, the 25 lh , and the 75'* percentiles, 
as a function of local density as defined above. We see an over- 
all shift of these distributions to lower values in high density 
regions. Furthermore, the most strongly star-forming galax- 
ies, i.e., the tail of the distribution for Ha EW > 5A, appears 
to experience the largest effect. To understand this effect fur- 
ther, we consider the (inverse) concentration index C of SDSS 
galaxies (Shimasaku et al. 2001). Although concentration will 



be sensitive to non-morphological features like the presence of 
nuclear star-formation, there appears to be a reasonable corre- 
lation between C and the Hubble morphological classifications 
as shown in Figure 10 of Shimasaku et al. (2001). For the anal- 
ysis in this paper, we have chosen a threshold of C > 0.4 to 
define late-type galaxies, instead of C > 0.33 as proposed by 
Shimasaku et al. (2001), as our choice provides a cleaner, but 
more incomplete, sample of such galaxies. We estimate from 
Figure 10 of Shimasaku et al. (2001) that our higher threshold 
ensures less than 5% contamination from early-type galaxies in 
the late-type sample. We find that the tail of the Ha EW dis- 
tribution is dominated by late-type (spiral) galaxies regardless 
of local density, i.e., > 75% of all galaxies with Ha EW > 5 A 
are consistent with being late-type galaxies. 

We note here that the 25'* percentiles of the Ha and [Oil] 
lines appear to follow different slopes with density, with the Ha 
25 th percentile showing almost no change with density. This is 
because the Ha EWs at < 5 A are dominated by stellar absorp- 
tion, and the strength of the absorption is only weakly depen- 
dent on the stellar population for galaxies older than a few Gyr. 
On the other hand, the [Oil] EW is unaffected by absorption 
and thus clearly shows the continual decline of emission flux 
into the densest regions. 

In Figure 4, we show the density-SFR relation for galax- 
ies that is analogous to the density-morphology relation of 
Dressier (1980) (although the density parameters are not di- 
rectly comparable). For this analysis, we have also removed 
obvious Active Galactic Nucleii (AGN) using the prescription 
outlined in Kewley et al. (2001) and using the [Nil], [Olll], Ha 
and H^ emission lines. If the detection of an AGN was am- 
biguous, based on the prescription of Kewley et al. (2001), 
then we do not exclude the galaxy from Figure 4. We also 
show the SFRN of galaxies as a function of local galaxy den- 
sity, which shows two interesting phenomena: First, when com- 
pared to the corrected SFR, there is a much larger difference be- 
tween the low and high density regions. This is likely an effect 
of the lower luminosity galaxies in the field having a higher 
rate of SFR per unit stellar luminosity than the more lumi- 
nous galaxies. Secondly, the "break", or characteristic density, 
seen in the density-SFR relation at a galaxy surface density of 
~ 1 h^ Mpc" 2 is much more prominent. At densities lower than 
this "break", the SFR of galaxies continues to increase, i.e., in 
Figure 4, the median and 75'* percentile of the SFRN distribu- 
tion keeps increasing all the way to the lowest densities studied 
in this paper (0.08 Mpc" 2 ). 

3.2. SFR as a Function of Clustercentric Radius 

We now turn our attention to virialized systems like clusters 
and groups of galaxies, and investigate the SFR of galaxies as 
a function of clustercentric radius. We test whether the fea- 
tures observed in the density-SFR relation, (e.g., the break or 
characteristic density at ~ 1 h^ 2 Mpc" 2 ) are also present in the 
radius-SFR relation around known virialized systems. The de- 
tails of how we objectively select clusters and groups from the 
SDSS EDR data are in the Appendix. We also discuss the 17 
clusters and groups used in this paper in the Appendix. 

In Figure 5, we present the observed distribution of Ha and 
[Oil] EW, as quantified by the median, 25'*, and 75'* per- 
centiles, as a function of (projected) clustercentric radius for 17 
groups and clusters of galaxies detected in the EDR data (see 
Appendix). To accommodate the wide range of masses spanned 
by our cluster sample, we have re-scaled all radial distances by 
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the virial radius (R v ) of each cluster. R v is computed from the 
line-of-sight velocity dispersion, in units of km s" 1 , using the 
formula R v ~ 0.002cr r ^ Mpc from Girardi et al. (1998). For 
this calculation, we used a r (l) as discussed in the Appendix. 
For consistency with the higher redshift work, we provide mea- 
surements for both the [Oil] and Ha emission lines. In Figure 
6, we show the same as in Figure 5, but now for the SFR and 
SFRN of galaxies as a function of (projected) clustercentric ra- 
dius. 

For comparison, we have also constructed a "non-cluster" 
(or field) sample of galaxies that consists of all galaxies within 
our volume-limited sample that are located, in redshift space, 
within 3.5<7 r of a cluster but are at a (projected) clustercentric 
distance of > 25 R v from any of our clusters. This methodology 
guards against potential redshift selection effects and/or red- 
shift evolution of the field population, because we have defined 
our field population in the same redshift shells as the cluster 
galaxies. We note however that the difference in the mean (or 
median) SFR of galaxies of the field population as defined us- 
ing our method and that of taking all galaxies regardless of the 
presence of clusters, is less than 10% (with our mean SFR mea- 
surement being systematically higher as expected). Our derived 
field values for the 25 th and 75 th percentiles are shown as lines 
on Figures 5 and 6. 

Figures 5 and 6 show a clear decrease in the star-formation 
activity (in either the SFR or the EW of Ha and [Oil]) as a func- 
tion of clustercentric radius. As seen in Figures 3 and 4, the 
SFR of galaxies in dense regions differ from that in the field in 
two ways: the whole distribution of SFRs (and EWs) is shifted 
to lower values, and the skewness of the distributions decreases, 
with the tail of the distribution containing high SFR galaxies di- 
minishing as one enters denser environments. As discussed in 
Section 3.1, this tail of strongly star-forming galaxies (Ha EW 
> 5 A) is dominated by late-type galaxies. For comparison with 
higher redshift studies, the mean (median) of the Ha and [Oil] 
EW distributions within one virial radius of the clusters studied 
herein are 3.4A (-0.7A) and 3.2A (1.8A), respectively. 

The results in Figures 5 and 6 are qualitatively similar to 
those for distant (z ~ 0.3) clusters of galaxies (Balogh et al. 
1997). A new aspect of our study, however, is that we can map 
this decrease in SFR all the way from the cluster cores into the 
field population. We have used the Kolmogorov-Smirnov (KS) 
statistic to test the distributions of EWs and SFRs, in each radial 
bin of Figures 5 and 6, against the distribution of EWs and SFRs 
derived from our field population. This test is designed to look 
for global differences between two distributions and allows us 
to determine the clustercentric radius at which the difference 
between the two distributions becomes statistically significant. 
The distribution of Ha EWs becomes different from the field at 
the 68% (> la) level at 3.2R V , while for [Oil] the distributions 
become different at 3.7^,. For SFR and SFRN, the effect of the 
cluster environment becomes noticeable, at the > la level, at 
clustercentric radii of 2.8R V and 2.4R V , respectively. 

It is evident from Figures 5 and 6 that the tails of these distri- 
butions are more strongly environment-dependent than the me- 
dians. Therefore, we have performed a separate statistical test 
on the tails of the distribution and how they change with clus- 
tercentric radius. We calculated the difference between the 25 th 
and 75 th percentiles in each radial bin of Figures 5 and 6, and 
the corresponding 25 lh and 75 th percentiles for the field (shown 
as solid, straight lines in Figures 5 and 6.) In order to assess 
the statistical significance of these differences, we created, for 



each radial bin, 1000 fake datasets via bootstrap re-sampling 
(with replacement) and computed, for each dataset, the same 
difference between cluster and field for both the 25 th and 75 th 
percentiles. This exercise provides a distribution of differences 
(between the cluster and field) for each radial bin for both the 
25 th and 75 th percentiles. We then determined the probability 
that the observed percentile difference is consistent with zero, 
i.e., we just count the fraction of percentile differences above 
zero in these fake datasets. Using this methodology, we deter- 
mined, at a greater than 68% confidence (> la), that the 75 th 
(25 th ) percentile of the Ha EW distribution becomes statisti- 
cally different from the field population at a clustercentric ra- 
dius of 3.5R V (0.5R t ). For the SFR distribution, the 75 th (25 th ) 
percentile becomes statistically different from the field popula- 
tion at a clustercentric radius of 3.6R V (0.5R V ). For the SFRN, 
the 75 th (25 th ) percentile becomes statistically different at a 
clustercentric radius of 3.9R V (0.9R^ 

These tests demonstrate that the SFR distribution of the 
galaxy population differs from the field (at la) out to ~ 3-4 
virial radii. These radii correspond to ~ Ah^\ Mpc, based on 
the average virial radii computed for the clusters given in Table 
1 . These tests confirm that all levels of SFR are affected, with 
the most strongly star-forming galaxies being affected more 
severely (i.e., the 75 th percentile) than the quiescent population. 
The implications of these results will be discussed in Section 4. 

To fully understand the correspondence between the radius- 
SFR relation and the density-SFR relation presented in Section 
3.1, we present the correlation between the local projected den- 
sity and clustercentric radius measured for each galaxy within 
the vicinity of our 17 clusters and groups. As expected, there is 
a linear correlation between these two quantities, and the char- 
acteristic density of ~ 1 Mpc -2 in Figure 4 corresponds to 
a range of clustercentric radii from ~ 2 to 3 virial radii. This 
range of radii is consistent with our observations above in that 
the distribution of SFRs in cluster galaxies begin to differ from 
the field population at ~ 3 to 4 virial radii. 

4. DISCUSSION 

4. 1 . Possible Biases and Systematic Errors 

Before we interpret our results, we must first investigate the 
effects of possible systematic biases on our work. These in- 
clude the effect of: i) The SDSS analysis software and the 
SDSS survey observing strategy; ii) Aperture bias due to the 
3 arcsecond fiber width; Reddening due to dust in the host 
galaxies; and iv) The smearing correction performed by the 
SDSS to improve the flux calibration of the continuum. We 
present a detailed analysis of these possible biases in the Ap- 
pendix, but we find that none of these potential biases have a 
large effect on our results. 

4.2. The Density-Morphology Relation 

As we move towards denser environments, the galaxy pop- 
ulation becomes dominated by early-type galaxies that have a 
lower intrinsic SFR (Kennicutt 1983; Jansen et al. 2000). In 
this section, we attempt to determine whether or not the SFR 
of galaxies of a given morphology are themselves affected by 
environment, i.e., is the density-SFR relation distinct from the 
density-morphology relation? Previously, Balogh et al. (1998) 
have attempted to answer this question using an objective mor- 
phological classification, based on a galaxy bulge-to-disk lu- 
minosity, that was then correlated with the radial dependence 
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of SFR in distant clusters. They concluded that cluster galax- 
ies of a given redshift, luminosity, and bulge-to-disk ratio had 
lower star formation rates compared with their counterparts in 
the field. Some evidence supporting this has also been reported 
by Hashimoto et al. (1998), Couch et al. (2001), Balogh et al. 
(2002a) and Pimbblet et al. (2001). We have attempted to ad- 
dress this issue in this paper in two different ways. First, using 
a subset of SDSS galaxies that possess both a visual morpho- 
logical classification (from either Dressier 1980 or Shimasaku 
et al. 2001) and a SDSS Ha EW measurement. Secondly, using 
objective morphologies based on the (inverse) concentration in- 
dex of SDSS galaxies as discussed in Section 3.1. We discuss 
these two tests below. 

We have cross-correlated our sample of galaxies with sam- 
ples of Dressier (1980) and Shimasaku et al. (2001). For the 
Dressier (1980) cluster sample, we find 114 galaxies in com- 
mon with our sample and thus have a SDSS Ha measurement. 
The sample includes 23 ellipticals, 39 spirals and 50 lenticulars 
that satisfy our selection criteria (Section 2.2). For the Shi- 
masaku et al. (2001) sample, which provides morphological 
classifications for 456 SDSS galaxies, selected over a represen- 
tative range of environments, we find 57 early-type galaxies 
and 73 late-type galaxies that satisfy our selection criteria (out- 
lined in Section 2.2). In this case, we have added together the 
ellipticals and lenticulars because of the concerns noted by Shi- 
masaku et al. (2001) regarding their tendency to preferentially 
classify lenticular as ellipticals, compared with the RC3 cata- 
log. Overall, these two samples of morphologically-classified 
galaxies cover a similar range in luminosity and redshift as the 
galaxy sample outlined in Section 2.2. The largest systematic 
uncertainty is in the consistency between the morphological 
classifications. 

In Figure 8, we show the distribution of Ha EWs for these 
galaxy samples as a function of their morphological type. We 
see a difference (confirmed using a KS test) in the Ha EW dis- 
tribution of spirals between these two samples, i.e., the Dressier 
sample of spirals has, on average, a lower SFR than seen in the 
Shimasaku et al. (2001) sample. This difference could be due 
to one or more of the following possibilities. First, it could re- 
flect a real environment effect, that cluster spirals have lower 
SFR than the average field spiral population. Second, it could 
be due to the fact that the morphological bins are too coarse to 
define a homogeneous galaxy population. A different relative 
distribution of Sa and Sc galaxies in the two samples, for ex- 
ample, could give rise to the same trend (see Shane & James 
2001). Third, the difference may reflect inconsistencies in the 
morphological classification criteria in the two samples. Thus, 
although these results are intriguing, a more definitive answer 
will require an automated classification of the SDSS galaxy 
morphologies, on a continuous scale. 

As a first attempt of this, we have used the (inverse) concen- 
tration index discussed in Section 3. 1 to investigate the degener- 
acy between the density-morphology and the density-SFR re- 
lations. This galaxy parameter can be used as an objective mor- 
phological classification, but is hard to interpret and to compare 
directly to the visual morphologies. In Figure 9, we present the 
distribution of SFRN as a function of projected local galaxy 
density (as originally presented in Figure 4), but split the sam- 
ple into two broad morphological bins based on the (inverse) 
concentration index, i.e., late-type galaxies (with C > 0.4, see 
Section 3.1) and early-type galaxies (with C < 0.4). As dis- 
cussed in Section 3.1, the tail of the SFR distribution (with 



high Ha EWs) is dominated by late-type galaxies at all den- 
sities. This plot also demonstrates that the late-type galaxies 
in our sample lie on a density-SFR relation similar to the one 
observed for the whole sample. As expected, the early-type 
galaxies have little, or no, star-formation, but they also appear 
to follow a shallow density-SFR relation (for this sample how- 
ever, we must be concerned about contamination from the late- 
type galaxies and the effects of stellar absorption). 

To quantify these effects in more detail, we present in Figure 
10, the distribution of Ha EW as a function of clustercentric ra- 
dius for all galaxies as well as just the late-type galaxies (with 
C > 0.4). We have performed a KS test between the late-type 
galaxies (C > 0.4) in the tails (> 5A) of the four Ha EW distri- 
butions in Figure 10 (the black filled histograms) and the late- 
type galaxies in the field (i.e., also with C > 0.4 and > 5 A). 
This analysis indicates that the inner three clustercentric radii 
(0.3R V , 0.9R V , 1 .6R V ) are inconsistent with the field population, 
at the 82.5%, 71% and 99.8% probability levels, respectively. 
This suggests that the distribution of Ha EWs for strongly star- 
forming (> 5 A) late-type galaxies does change with density, 
which is consistent with our results in Section 3.1 and with 
those of Hashimoto et al. (1998). 

4.3. Physical Interpretation 

The most striking result of this paper is the critical density 
(or radius) where the SFR of galaxies changes from that of the 
field. This happens at a density of ~ 1 h^l Mpc" 2 , which cor- 
responds to a clustercentric distance of between 2 and 3 virial 
radii. In this section, we attempt to understand this observation 
within the context of hierarchical structure formation. In hier- 
archical models, clusters and groups form from slightly over- 
dense regions in the initial matter distribution. At first, clusters 
expand with the Hubble flow until they reach a maximum size 
and then de-couple from the expansion and collapse to form a 
virialized object. Using the spherical collapse model (Gunn & 
Gott 1972), one can show that the point of maximal size, also 
called the "turn-around" radius (R t ), is equal to a density of 
approximately six times the mean density of the background at 
that time. The virial radius (R^ corresponds to a radius where 
the density is a few hundred times the mean density of the back- 
ground. Therefore, there are two natural scales for an overden- 
sity in the Universe - the virial radius and the turn-around ra- 
dius - and the latter radius marks the limit of the gravitational 
influence of a cluster or group. The region between the two 
scales is known as the infall region. 

We can estimate the ratio of virial (R v ) and turn-around (R,) 
radius as follows. Suppose that the region between R v and R t 
is completely empty. Then M t jR] = M v /Rl (R v /R,) 3 , where M, 
and M v are the masses of the cluster within the turn-around and 
virial radii, respectively. Using M t jR] ~ 6 and M v jR\ ~ 300, 
we find R, ~ 4R V . In practice however, the region between 
R v and R, is not empty and a more realistic case is to as- 
sume that clusters or groups of galaxies follow, for example, 
the NFW mass density profile (Navarro et al. 1996). We can 
thus find the turn-around radius by integrating the profile until 
M(<R)/R 3 = 6. In this case, we find that R t ~ 5 R v , which is in 
good agreement with our naive calculation above. 

It is intriguing that this scale (R, ~ 5R V ) is comparable to the 
scale at which the cluster SFR becomes statistically indistin- 
guishable from the field. However, it is unclear that this radius 
has any physical meaning that might effect the SFR of indi- 
vidual galaxies, as the mechanisms proposed for changing the 
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SFR of galaxies in the cores of rich clusters, e.g., ram-pressure 
stripping of the gas, galaxy harassment (Moore et al. 1999), 
tidal disruption (Byrd & Valtonen 1990), are unlikely to be im- 
portant at low densities and/or large virial radii. A more ap- 
propriate mechanism for changing the SFR of these galaxies 
may be the merger or close tidal interactions of galaxies in less 
dense groups within the infall regions of clusters (Zabludoff 
et al. 1996; Zabludoff & Mulchaey 1998; Kodama et al. 2001). 

To obtain a more detailed physical understanding of our re- 
sults, we must compare our observations with numerical simu- 
lations. For example, Balogh et al. (2000) and Diaferio et al. 
(2001) have recently showed that simple (heuristic) models 
for the gas depletion of galaxies within hierarchical models of 
structure formation are able to accurately reproduce the de- 
crease in the star-formation rate of galaxies seen within the 
cores of CNOC1 clusters of galaxies (Balogh et al. 1997, 1998). 
In both these works, the SFR of a galaxy is affected by both lo- 
cal and external processes. Locally, the SFR is governed solely 
by the consumption rate of cold gas in the disk, dependent only 
on gas density and the feedback model. The only external pro- 
cess that effects the SFR of a single galaxy is the stripping of 
its hot gas reservoir after it merges with a larger (e.g., group 
or cluster) halo. Following this, the SFR declines gradually, 
as the galaxy consumes the remaining cold, disk gas. There- 
fore, the main physical properties of a galaxy in the simulation 
that controls these processes are the amount of cool gas in the 
galaxy and the time since the last interaction with a larger halo. 
Diaferio et al. (2001) predicts that the mean SFR of galaxies in 
clusters should be lower than the field out to 2R V (here R v is a 3- 
dimensional radius from the simulations, while throughout the 
paper, we have quoted R v as a projected radius from the cluster 
cores). This reduction is due to galaxies that had been near the 
cluster core, but thrown out to large radii during major merg- 
ers (Balogh et al. 2000; Evrard & Gioia 2002). Our results are 
qualitatively similar to the predictions of Balogh et al. (2000) 
and Diaferio et al. (2001) in that these hierarchical models of 
structure formation can affect the SFR of galaxies beyond the 
virial radius. However, we will require further comparison of 
the simulations and observations to determine if these simple 
models are all that is required to explain the data, or whether 
we need additional physical processes that effect the SFR of 
galaxies at larger radii in the infall regions of clusters, i.e., from 
2R V to 4R V . 

4.4. Comparison with Previous Work 

A similar analysis has been carried out using the LCRS 
(Hashimoto et al. 1998) and 2dF Galaxy Redshift Survey (2dF- 
GRS) data (Lewis et al. 2002). For the 2dFGRS, the luminosity 
limits, redshift limits, and Ha EW distributions are comparable 
with our sample, and both studies are remarkably consistent in 
their main conclusions. First, star formation is reduced, rela- 
tive to the field, out to > 3 virial radii from known clusters and 
groups. Secondly, there appears to be a critical density of 1 
galaxy (brighter than M* + 1 .5) per Mpc 2 , below which there is 
a weaker correlation with star-formation rate. 

Our work moves beyond the 2dF study in several ways. Most 
importantly, we have computed a local density for every galaxy 
in the EDR, regardless of its proximity to a cluster. More- 
over, our cluster catalog is dominated by systems with much 
lower velocity dispersions than those in Lewis et al. (2002), so 
it further emphasizes that the density-SFR relation is universal 
and does not depend on the large-scale mass of the embedding 



structure. Furthermore, Lewis et al. (2002) are unable to com- 
pute Ha fluxes, and hence star formation rates, due to uncer- 
tainty in the continuum fluxing. This is not a problem for the 
SDSS, and we have been able to show that there exists an ab- 
solute density-SFR relation, as well as a relative one. We have 
also used the fluxes from another line ([Oil]) to demonstrate 
that our results are robust to the effects of reddening. Finally, 
the data analysed by Lewis et al. (2002) are preliminary in the 
sense that sampling is not complete over their full extracted re- 
gions. The consistency between our results shows that this is 
not a large problem. 

Recently, Kodama et al. (2001) reported the detection of a 
"break" in the colors of galaxies around distant cluster Abell 
851 (z = 0.41). They observed an abrupt change in the col- 
ors of galaxies at a local surface density of galaxies of ~ 
20011^5 Mpc" 2 , which corresponds to a radial distance of ~ 
1 h^ Mpc from the center of this cluster. Even after correcting 
for the differences in the luminosity limits of the two data sets 
(0.025Ly in Kodama et al. (2001) compared to 0.4L*« here), 
the "break" seen by Kodama et al. (2001) occurs at an order of 
magnitude greater surface density, and also a smaller cluster- 
centric radius, than the "break" reported here. Further observa- 
tions will be needed to determine if this difference is physical 
(i.e., due to differences in the redshift or luminosity limits of 
the two studies) or an artifact of the different analyses. 

4.5. Future Work 

As the SDSS database increases, we will be able to increase 
the size and depth of the sample used herein, i.e., to increase 
the number of galaxy clusters, as well as to relax some of the 
conservative selection criteria we have imposed on our present 
sample. Furthermore, we will use additional information on 
these galaxies extracted from both the SDSS photometry and 
spectra. For example, we can use the five passband SDSS pho- 
tometry to look for color gradients within these galaxies to fully 
correct for the aperture biases as well as study how the envi- 
ronment effects the star formation within a galaxy (see Moss 
& Whittle 1993, 2000; Rose et al. 2001). We must also ex- 
tend our study to characterize the effects of the environment on 
the overall star formation history of galaxies. For instance, we 
should include post-starburst galaxies (E+A, k+a and a+k; see 
Zabludoff et al. 1996; Poggianti et al. 1999; Balogh et al. 1999; 
Castander et al. 2001), because these may be key to fully un- 
derstanding star-formation activity within dense environments. 
In addition, we should probe fainter luminosities to determine 
the rate of change of SFR per unit stellar luminosity. 

Most importantly, the research discussed in Section 4.2 
clearly requires a more robust, automated morphology for each 
of the SDSS galaxies similar to the bulge-to-disk decomposi- 
tions used by others (Balogh et al. 1998; Balogh et al. 2002b). 
We can then perform a multi-variate analysis on the galaxy 
population that treats the SFR, local density, morphology, mass 
and luminosity of the galaxies in a self-consistent manner. We 
must also expand our definition of local density to include more 
information on the higher-order moments of the galaxy distri- 
bution, e.g., filaments. A full characterization of the effect of 
the environment on the overall star formation history of galax- 
ies is paramount to our understanding of the physical processes 
that could be responsible for galaxy evolution. 
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5. CONCLUSIONS 

We have investigated the effects of the local galaxy environ- 
ment on the star-formation of galaxies using the Early Data 
Release of the Sloan Digital Sky Survey. For this work, we 
have restricted our analysis to a volume-limited sample of 
8598 galaxies brighter than M(r*) < -20.45 (k-corrected, for 
H = 75kms _1 Mpc" 1 ) over the redshift range 0.05 < z < 0.095. 
For all galaxies, we have characterized the star-formation us- 
ing the Ha and [Oil] EWs, as well as by computing their 
star-formation rates (SFR) and normalized star-formation rates 
(SFRN; SFR per unit luminosity) from the SDSS data. We have 
quantified the local galaxy environment using the projected dis- 
tance to the 10 th nearest neighbor brighter than M*+ 1 as well 
as by measuring the clustercentric distance of galaxies from the 
cores of known clusters and groups of galaxies (see Appendix). 
We have extensively tested our results as discussed in detail in 
the Appendix. 

This work expands upon previous studies of the SFR of 
galaxies as a function of environment in four ways: i) We ex- 
tend such studies into the group and poor cluster regime; ii) 
We focus on the low redshift universe; Hi) We have traced the 
SFR activity of galaxies well beyond the cores of the clusters 
and groups into the field population; and iv) We are able to ac- 
curately quantify the local galaxy density in a uniform manner 
that is not subject to statistical background correction. We con- 
clude that: 

• The distributions of [Oil] and Ha EWs, as quantified by 
the median and 25 th and 75"' percentiles, changes as a 
function of local galaxy density as measured using the 
distance to the 10 th nearest neighbor. We witness simi- 
lar relations in the SFR and SFRN of galaxies as shown 
in Figure 4. This effect is characterized in three ways. 
First, there is shift in the overall distributions of EW, 
SFR, and SFRN to lower values with increasing galaxy 
density. Secondly, the skewness of the distributions de- 
creases with increasing local galaxy density, i.e., the tail 
of strongly star-forming galaxies (Ha EW > 5A), as 
quantified by the 75'* percentile of the distribution, is 
noticeably decreased in high density regions. Finally, 
we see a "break" (or characteristic scale) in the corre- 
lation between SFR (and SFRN) and density at a local 
galaxy density of ~ 1 Mpc" 2 (for galaxies brighter 
than M*). Figures 3 and 4 represent the density-SFR 
relation of galaxies that is analogous to the density- 
morphology relation of galaxies (Dressier 1980). 

• The distribution of EW, SFR, and SFRN of galaxies, as 
quantified by the median and 25 th and 15 th percentiles, 
changes as a function of clustercentric radius for the 17 
clusters and groups discussed in this paper (see Table 
1). This effect is most noticeable for the strongly star- 
forming galaxies in the 15 th percentiles of the distribu- 



tion. Using a KS test and boot-strapping re-sampling 
techniques, we find that the star formation rates of 
galaxies begin to decrease, compared with field galax- 
ies, starting at 3-4 virial radii (with > ler statistical sig- 
nificance). Within one virial radius of our clusters, the 
means (medians) of the Ha and [Oil] EW distributions 
are 3.4A (-0.7A) and 3.2A (1.8A), respectively. 

• As shown in Figure 7, the break or characteristic density 
seen at a (projected) galaxy density of ~ 1 Mpc in 
the density-SFR relation corresponds to ~ 2 and 3 virial 
radii for the systems discussed in this paper. Therefore, 
there is good agreement between our results from the 
clusters and groups (discussed in Section 3.2) and the 
general density-SFR relation (discussed in Section 3.1). 

• We have investigated the possible degeneracy between 
our density-SFR relation and the density-morphology 
relation (Dressier 1980). There is some evidence that 
the SFR of galaxies in dense regions varies even among 
galaxies of the same morphology. For example, we have 
used the (inverse) concentration index of SDSS galax- 
ies to show that the tail of the strongly star-forming 
(Ha EW > 5 A) galaxies observed in the Ha EW, SFR, 
and SFRN distributions is dominated (> 75%) by late- 
type (spiral) galaxies. Using a KS test, we have shown 
that this tail of late-type galaxies is different in dense 
regions (within 2 virial radii) compared with similar 
galaxies in the field. We stress, however, that the sig- 
nificance of these results remains uncertain at present 
due to potential systematic biases in the morphological 
classifications used herein. 

• Our results are qualitatively consistent with the predic- 
tions of Balogh et al. (2000) and Diaferio et al. (2001) 
in that these hierarchical models of structure formation 
can affect the SFR of galaxies well beyond the virial ra- 
dius. However, we require more detailed comparisons 
between the simulations and observations to determine 
if such simple (heuristic) models are all that is neces- 
sary to explain the data, or whether we need additional 
processes that affect the SFR of galaxies at large clus- 
tercentric radii and/or low densities. 

• Our results are in good agreement with the recent 2dF 
work of Lewis et al. (2002) as well as consistent with 
previous observations of a decrease in the SFR of galax- 
ies in the cores of distant clusters (Balogh et al. 1997; 
1998). Taken together, these results demonstrate that the 
decrease in SFR of galaxies in dense environments is 
a universal phenomenon over a wide range in densities 
(from the rarefied field to poor groups to rich clusters) 
and redshifts. 
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APPENDIX 

Selection of the Groups and Clusters of Galaxies 

A critical part of our study of galaxy star-formation activity in dense environments is the selection of virialized groups and clusters 
of galaxies within the EDR. This was achieved using the objective C4 algorithm described in Nichol et al. (2001) and Miller et al. 
(2002, in prep.), which assumes that galaxy groups and clusters have a co-evolving set of galaxies with similar colors, e.g., within 
the E/S0 ridgeline. Like any cluster-finding algorithm, this does bias our selection of systems because we assume a model for the 
objects that we are searching for. However, such biases can be quantified using simulations (see Miller et al. 2002, in prep.) and 
is preferred to using visually-compiled catalogs, like Abell (Abell, Corwin, & Olowin 1989), whose selection function is hard 
to quantify through simulations. We stress that these criteria are used only to find the clusters and not to assign membership for 
subsequent analysis. 

We briefly outline the methodology we used herein which is based on the C4 algorithm discussed in more detail in Miller et al. 
(2002, in prep). For each galaxy in our sample, we count the number of neighbors within a seven-dimensional box (four colors, RA, 
DEC and a redshift). The size of the box in the color dimensions is determined from the errors in the galaxy magnitudes and the 
intrinsic scatter in galaxy colors within the cores of groups and clusters. The angular size corresponds to 1 h^Mpc at the redshift of 
the galaxy. The size of the box in the redshift direction is Az = 0. 1 and thus only rejects obvious foreground and background galaxies. 
We then compare this galaxy count (around each galaxy) to the distribution of neighbor counts for the whole field population and 
calculate the probability that the target galaxy is a field galaxy. The field distribution is determined by moving the seven-dimensional 
box to 100 randomly chosen galaxies with similar seeing and galactic extinction measurements. Using the False Discovery Rate 
thresholding technique of Miller et al. (2001), we then identify galaxies that have a low probability of being in the field, i.e., such 
that < 5% of these galaxies could be mistaken field galaxies. We note that the C4 algorithm does not assume a given color for 
these galaxies; it only assumes that these galaxies have similar colors. This non-parametric approach significantly reduces the risk of 
biasing our cluster sample because we are not modeling the expected colors of the cluster galaxy population (Gladders & Yee 2000; 
Goto et al. 2002). The details of the full C4 algorithm, including the selection function and completeness limits, are presented in 
Miller et al. (2002, in prep.). 

We first applied the above implementation of the C4 methodology to the EDR spectroscopic data, i.e., to all galaxies with a redshift. 
This finds ~ 5% of galaxies to be located in highly clustered regions (in our seven-dimensional space), and it is therefore possible to 
locate candidate clusters and groups within this highly-clustered galaxy dataset using a non-parametric density estimator. For each 
candidate cluster, we then compute the mean centroid and redshift of each system using all of the EDR data available, i.e., we now 
include all galaxies in the EDR spectroscopic and photometric data sets regardless of the C4 algorithm. In addition, we also estimate 
the line-of-sight velocity dispersion (o>) of each system using a robust o> estimator (see Miller et al. 2002 in prep.) that iteratively 
rejects all galaxies > 3o> away from the mean redshift. We then fit a Gaussian to the velocity distribution of each candidate cluster 
and keep only systems for which the difference between the standard deviation about the mean and the la dispersion of the fitted 
Gaussian is less than a factor of two. This approach removes spurious systems that account for less than 10% of whole sample and 
does not affect our results. 

We present in Table 1 the 17 groups and clusters of galaxies that were found in the EDR via this implementation of the C4 
algorithm and that satisfy the selection criteria discussed above (in total, 30 clusters were detected in the EDR, but 13 of these were 
outside the redshift limits used herein). We provide the RA and DEC of the cluster centroids (in degrees), their mean redshift, the 
line-of-sight velocity dispersion derived from the Gaussian fit to the spectroscopic data (o>(l)), the line-of-sight velocity dispersion 
derived from the standard deviation of the galaxy distribution (o>(2)), the number of cluster members within 2 virial radii (h = 0.75) 
of the cluster center, brighter than M r = -20.45, and within ±3a>(l) of the mean cluster redshift, and the X-ray luminosity derived 
from the ROSAT All-Sky Survey (RASS, Voges et al. 2000) data in units of lO^ergs" 1 (0.5-2.4 keV; h = 0.75). The last column 
of the table provides the name (and redshift if available) of the nearest known cluster within 15 arcminutes from the NED database. 
Details of these systems are discussed in Miller et al. (2002, in prep.) along with redshift histograms and color-magnitude relations 
that illustrate the reality of these overdensities. 

As can be seen from Table 1, our C4 systems span a wide range of velocity dispersions (from c± 200kms _1 to 1000 kms" 1 ) and 
X-ray luminosities (from lO^/z^ergs" 1 to <~ 4 x lO^/^fergs -1 ). The relation between these two physical quantities is consistent 
with the L x -o> relation of Zabludoff & Mulchaey (1998) and Mahdavi & Geller (2001), while the ratio of groups to clusters in our 
sample is consistent with expectations based on the volume of our galaxy sample, i.e., in c± 460 deg 2 of sky, we would only expect to 
find a few X-ray luminous clusters (Bohringer et al. 2001). Therefore, our sample mostly covers the group and poor cluster regime 
(~ 10 43 hjjerg s" 1 in X-ray luminosity). This is different from most previous studies of environment-dependent star-formation which 
have focused on rich clusters, and thus extends such studies to more common galactic environments. 

SDSS Analysis Pipelines, Survey Strategy and Design 

As discussed in Stoughton et al. (2002) and Frieman et al. (in prep.), the SDSS data analysis pipelines have been extensively 
tested and can accurately determine the main characteristics of emission and absorption lines, especially at the signal-to-noise ratios 
used herein. As part of a detailed check on these measurements, we have studied the 1231 duplicate observations of SDSS galaxies 
discussed in Section 2.2. In Figure 11, we show the percentage difference between these duplicate observations of Ha EW as a 
function of EW. As expected, the measurement error is large for small EWs, but approaches ~ 10% for Ha EW> 5 A. We further 
propagate these errors through to our measurements of the SFR, and present in Table 3 the mean and median percentage difference 
of the SFR as a function of SFR. For low SFRs, the error is substantial, but above a SFR of ~ 1 M Q yr -1 , the mean and median error 
is below 20%. 
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Stellar Absorption 

The present SPECTROID analysis pipeline does not account for possible stellar absorption, and we have made no correction 
for this in our measurements of the Ha emission line fluxes in this paper. This does not affect our results because: i) The effect 
is expected to be < 5 A, which is small relative to the strong emission-line galaxies in which we are interested; ii) We include all 
galaxies (regardless of the size and sign of their measured EW) when calculating the percentiles of the EW and SFR distributions. 
If stellar absorption is only a weak function of EW or SFR, then its effect will be to systematically shift the whole distribution to 
smaller values that, in turn, will be reflected in the measured median, 25'*, and 75 th percentile of the EW distribution; Hi) Our results 
are consistent for both the [Oil] and Ha emission lines, which is reassuring because the [Oil] emission lines are unaffected by stellar 
absorption. Furthermore, [Oil] and Ha reside in different parts of the SDSS spectrum (and were observed on different CCDs in the 
SDSS spectrographs) and thus suffer different potential problems, e.g., the determination of the continuum near [Oil] can be severely 
affected by the 4000A break. Therefore, the fact these two lines give similar answers suggests that our results are not significantly 
affected by stellar absorption. 

Aperture Bias 

Aperture bias is a concern because we have observed all our galaxies through a fixed angular (3 arcseconds) fiber that is smaller 
than the angular size of our galaxies (see Konchanek, Pahre & Falco 2000 for a discussion of possible effects of aperture bias). 
Such a bias could result in a systematic increase in the observed Ha EW and SFR for higher redshift galaxies relative to lower 
redshift galaxies, because more of the galaxy light is passing down the fiber. To minimize this potential bias, we have followed the 
conclusions of Zaritsky, Zabludoff & Willick (1995) and restricted our sample of galaxies to z > 0.05. Above this redshift, Zaritsky 
et al. (1995) showed that for the LCRS, the spectral classifications of LCRS galaxies, which were observed through a 3.5 arcsecond 
fiber, are statistically unaffected by aperture bias. In fact, the results of Zaritsky et al. (1995) suggest that only ~ 10% of spiral 
galaxies in our sample are mis-classified as non-star-forming due to the inability of the fiber to capture Ha emission in the outer 
disk of the galaxy. 

To study aperture bias in our SDSS galaxy sample, we have split the sample discussed in Section 2.2 into three subsamples as 
a function of redshift: 0.05 < z < 0.072, 0.072 < z < 0.085, and 0.085 < z < 0.095. We chose these redshift shells because they 
represent equal volumes. In Table 2, we show the median and 75'* percentile of the Ha EW distribution for these three redshift shells 
as a function of galaxy type. First, we analysed all galaxies in each shell and find no evidence for a systematic increase in the median 
and 75 ,h percentile with redshift. Next, we restricted the analysis to just late-type galaxies using the (inverse) concentration index 
because we expect such galaxies to have the largest aperture bias. Once again, we find that the median and 75"' percentile of the Ha 
EW distribution of these galaxies does not significantly increase with redshift. As a final check for aperture bias, we further split the 
sample of late-type galaxies discussed above using their Petrosian g* - r* color and physical size (as defined by the Petrosian radius 
that contains 90% of the galaxy light in the r* band), and search for the largest observed increase in the median and 75 th percentile of 
the Ha EW distribution with redshift. The result of this search is given in Table 2 (labelled "Worst Case") where we find the smallest 
(< 8 hjl), bluest (g*-r* < 0.75), late-type galaxies possess the largest systematic increase in the median and 75 ,h percentile of the 
Ha EW distribution with redshift. This is reasonable because for these intrinsically small late-type galaxies, a larger fraction of the 
light captured by the fiber is from the disk. 

Thus the magnitude of any aperture bias is, at worst, 6A in the Ha EW and only affects a small fraction of the late-type galaxies 
(see Table 2). We use a KS test to compare the distribution of local (projected) densities observed in the three redshift shells for the 
subsample of small, blue, late-type galaxies discussed above, and find no statistical evidence for any difference in these distributions. 
This result indicates that this subsample of galaxies inhabits the same environments at all redshifts and that our conclusions are not 
significantly affected by any aperture bias. 

We note here that the recent work of Moss & Whittle (2000) and Bartholomew et al. (2001) both suggest that star-formation in 
cluster galaxies may be more centrally concentrated than in field galaxies. If true, this could introduce an aperture bias but in the 
wrong direction to explain our results, i.e., it would reduce the contrast between the cluster and field SFR distributions, because we 
would be underestimating the SFR in field galaxies by a larger factor than for cluster galaxies. Therefore, the density-SFR trends we 
observe will only be enhanced if this effect is important. 

Reddening 

Galactic extinction has been corrected for using the Schlegel et al. (1998) maps as discussed in Stoughton et al. (2002), while no 
correction has been made for possible intra-cluster extinction because this is expected to be small or non-existent (Bernstein et al. 
1995; Stickel et al. 2002). More worrisome, however, is internal dust extinction in each of the sample galaxies and its effect on the 
observed SFR (Hopkins et al. 2001; Sullivan et al. 2001; Chariot et al. 2001). To address this possible bias, we present in Figure 12 
the median and 75 th percentile of the SFR as a function of density, for three different SFR prescriptions. The first is the "standard" 
derivation, based on an empirical methodology that uses the Kennicutt (1998b) relation as given in Eqn 1., assuming one magnitude 
of extinction for all galaxies. The second, which is the prescription we have adopted throughout this paper, is the same but corrected 
for SFR-dependent reddening using the empirically derived relations of Hopkins et al. (2001). Finally, we use the direct modeling 
approach of Chariot & Longhetti (2001), who provide a more complete model for the dust, gas and ionizing radiation in galaxies 
based on a combination of emission lines. However, the Chariot & Longhetti prescription requires that the [Oil] , Ha , Hp, [Olll], 
[Nil] and [Sll] rest-frame EWs are larger than 5 A, which demands high signal-to-noise data and the presence of strong emission- 
lines. As a result, only 85 galaxies satisfy the Chariot & Longhetti (2001) prescription. For the remainder we just use the corrected 
SFR of Hopkins et al. (2001). Therefore, in Figure 12, the corrected SFR and Chariot & Longhetti (2001) SFR are very similar. 
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The fact that all three models show a trend of SFR with density suggests that the dust properties of galaxies do not vary sufficiently 
with radius to give rise to the observed trends. Finally, in Figure 3, we show that the [Oil] EW distribution as a function of density is 
similar to the result obtained using Ha . This is encouraging because the Ha and [Oil] lines have very different dust and metallicity 
sensitivities, come from different star-forming regions in the galaxies, and are separated by nearly 3000A in wavelength. 

Flux Calibration 

As discussed in Stoughton et al. (2002), a 4-minute smear exposure is taken for every SDSS spectroscopic science observation. 
This short exposure is designed to obtain a more accurate spectrophotometric calibration as the telescope is moved slightly during the 
smear exposure, effectively increasing the nominal 3 arcsecond fiber aperture to a 5 by 8 arcsecond 2 aperture. This produces a low 
signal-to-noise spectrum of the galaxy, but for a larger aperture, that then allows us to better correct the spectral continuum radiation 
for the affects of seeing and atmospheric refraction. Emission lines are excluded from the smear spectrum, so only the continuum 
of the science exposure is effected by this procedure. This could introduce a systematic bias into the EW measurements because the 
smear correction is only changing the continuum shape relative to the unchanged emission line properties. Smear exposures are now 
a standard part of the SDSS observing strategy. However, for the EDR, only half of the data were taken with a smearing exposure 
(see Table 3 in Stoughton et al. 2002) so we have a unique opportunity to study the effect of this observing strategy on the line 
measurements of SDSS galaxies. Using a Kolmogorov-Smirnov test, we were unable to detect any statistically-significant difference 
between the two distributions (smeared and non-smeared data) for Ha EW, [Oil] EW and the SFR. This test suggests that the smear 
correction does not systematically bias these measurements. 

Monte Carlo Simulations 

As a final test of our results, we performed 100 bootstrap simulations of our data which involved randomly re-assigning (with 
replacement) the SFR of each galaxy in our sample. The results of this test are consistent with no observed correlation between 
SFR and local galaxy density. Thus, the observed correlation between SFR and density in our data is not an artifact of the geometry 
or design of the sample because such systematic effects (minimum separation of fibers, boundary effects) are still present in our 
bootstrap simulations and would thus still appear as a correlation between SFR and density if significant. 
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Table 1 



Groups and clusters of galaxies used in this paper. 



RA (J2000) 


DEC (J2000) 


z 


ov(l) 


<J r (2) 


N 


L Y (44)h;l 


Name 








km s — ^ 


km s — ^ 








5.853 


-0.857 


0.063 


225 


430 


17 


0.02 




7.304 


-0.199 


0.059 


345 


515 


24 


0.16 


ZwCl 0027.0-0036 


15.574 


-0.379 


0.050 


312 


246 


7 


0.07 




20.039 


0.032 


0.077 


223 


351 


15 


0.04 




22.852 


0.627 


0.079 


504 


761 


39 


0.39 


A208 (z=0.0798) 


23.636 


-0.675 


0.079 


409 


592 


42 


0.16 




24.424 


-0.440 


0.056 


298 


556 


18 




ZwCl 0134.8-0045 


30.990 


0.270 


0.076 


228 


254 


10 


0.04 


A299 


151.942 


0.496 


0.096 


542 


563 


64 


1.23 


A933 (0.0956) 


191.831 


0.143 


0.088 


1243 


1450 


236 


1.96 




197.967 


-0.751 


0.084 


878 


1003 


85 


3.75 




199.787 


-0.873 


0.083 


677 


771 


67 




ZwCl 1316.4-0044 


201.575 


0.232 


0.082 


468 


528 


36 


1.07 


RXJ1326.3+0013 
(z=0.0820) 


206.317 


0.098 


0.088 


345 


330 


32 




ZwCl 1342.4+0016 


217.565 


0.337 


0.055 


297 


354 


14 






227.302 


-0.250 


0.090 


682 


701 


84 


0.34 


A2026 (z=0.0876) 


227.734 


-0.103 


0.091 


567 


573 


57 




A2030 (z=0.0919) 



We have included in our analysis clusters marked with a * even though their mean redshift is beyond the nominal redshift limit of the sample defined in the text. This 

is because some of their cluster galaxies are in the sample limits. 



Table 2 

Median and 75" ! Percentile of the Ha EW Distribution as a Function of Redshift and Galaxy Type 





0.05 < z < 0.072 


0.072 < z < 0.085 


0.085 < z < 0.095 




Median, 75 th (Numb.) 


Median, 75' h (Numb.) 


Median, 75'* (Numb.) 


All Galaxies 


8.1, 12.3(1632) 


7.9, 12.8 (2849) 


8.4, 13.1 (1977) 


Late-type 


15.3,20.4 (632) 


16.1,23.3 (1045) 


16.0, 23.1 (809) 


Worst Case 


27.1,35.8 (130) 


28.5,36.1 (205) 


31.5,41.8 (133) 



Table 3 

The Mean and Median Percentage Difference in SFR for Duplicate Observations as a function of SFR. 



SFR ( Mq yr -1 ) 


Median 


Mean 


<0.05 


25% 


141% 


0.05-0.15 


36% 


64% 


0.15-0.90 


16% 


29% 


> 3.0 


18% 


19% 
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FIG. 1 . — (Top) The distribution of observed Ha EW for the galaxies in our sample. The star symbol is the median value, while the triangle symbol is the mean. 
The negative EWs in the distribution are caused by statistical uncertainties and stellar absorption (see text for discussion). (Bottom) The distribution of observed 
[Oil] EW for the galaxies in our sample. The star symbol is the median value, while the triangle symbol is the mean and the diamond represents the mean field [Oil] 
EW from Balogh et al. (1997). The negative EWs in the [Oil] EW distribution are caused by statistical scatter. 
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FIG. 2. — The distribution of local galaxy densities (A 7 | Mpc 2 ) for galaxies in our sample. 
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FIG. 3. — (Top) The shaded area represents the distribution of Ha EW as a function of the projected local surface density of galaxies. The top of the shaded area 
is the 75 ,h percentile of the Ha EW distribution, while the bottom is the 25 < h . The median is shown as a line. We have used all available galaxies in the SDSS EDR 
that satisfy our selection criteria regardless of their location near a known cluster of galaxies. We have, however, excluded galaxies close to the survey boundary 
(see the text for a complete explanation) and those that may have an AGN (see text). Each bin contains 250 galaxies. (Bottom) The same as above, but for the [Oil] 
emission line. Each bin contains 150 galaxies. 
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FIG. 4. — (Left) The shaded area represents the distribution of corrected SFR (Hopkins et al. 2001) as a function of the projected local surface density of galaxies. 
(Right) The shaded area represents the distribution of SFRN (the normalized SFR; see text) as a function of the projected local surface density of galaxies. In both 
plots, the top of the shaded area is the 75 ,h percentile, while the bottom is the 25 ,h percentile. The median is shown as a solid line. We have used all available 
galaxies in the SDSS EDR that satisfy our selection criteria. We have excluded galaxies near the edge of the survey and those which may have an AGN present 
based on the Kewley et al. (2001) prescription. Each bin contains 150 galaxies. These plots represents the density-SFR relation that is analogous to the density- 
morphology relation of Dressier (1980). We note here that the SFRs presented here are not corrected for the 3 arcsecond SDSS fiber aperture and are therefore, 
systematically lower, by a factor of ~ 5, compared to total SFRs derived from the radio or by integrating the light from the whole galaxy (see Hopkins et al., in prep) 
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FIG. 5. — (Left) The shaded area is the distribution of Ha EW as a function of (projected) clustercentric radius. We have scaled the projected radial distances by 
the virial radii for the appropriate cluster ay(l) (see Table 1). The top and bottom of the shaded area represents the 75 '* and 25 ,h percentile of the EW distribution, 
while the light grey line inside the shaded area is the median of the distribution. The straight lines at the top and bottom of the shaded area show the 75 th and 
25 ,h percentile of the EW distribution for the field population respectively (see text for definition of the field population). (Right) The same, but for the [Oil] EW 
distributions. In each plot, there are 180 galaxies per bin. 
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FIG. 7. — The clustercentric distance as a function of the local galaxy density for galaxies in our cluster sample. The points are the median of each bin with the 
error on the median from bootstrap re-sampling. 
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FIG . 8 . — (Top) We show the observed distribution of Ha EWs as a function of morphological classification for 1 14 cluster galaxies serendipitously observed by 
the SDSS. The morphologies come from Dressier (1980) and include ellipticals (E), lenticulars (SO) and spirals (S). The Ha EWs come from the SDSS database. 
(Bottom) The distribution of Ha EWs for the Shimasaku et al. field galaxy sample. We only show late (S) and early-type (E/S0) galaxies (see text). 



22 



Gomez et al. 




FIG. 9. — (Left panel) The same as shown in Figure 4. (Middle panel) The density-SFR for late-type galaxies (i.e., with C > 0.4 as discussed in Section 3.1). As 
in other plots, the top of the shaded area is the 75 ,h percentile of the distribution, while the bottom of the shaded region is the 25 ,h percentile of the distribution. 
The median is shown as a black line near the middle of the shaded region. (Right Panel) The density-SFR for early-type galaxies (i.e., with C < 0.4). 
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FIG. 10. — The distribution of Ha EWs as a function of clustercentric radius and galaxy (inverse) concentration index. The unfilled histograms represent the whole 
distribution of Ha EWs regardless of morphological classification. The filled histograms are the distribution of Ha EWs for galaxies with an (inverse) concentration 
index of C > 0.4 (Shimasaku et al. 2001) and are thus classified as late-type galaxies. 
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FIG. 1 1 . — We show here the percentage absolute difference for the EW of Ha as function of the observed Ha EW, for EDR galaxies that were observed twice 
in the SDSS and analysed independently through the official SDSS SPECTROID analysis pipeline. The dashed vertical line is at Ha=5A, where the errors begin to 
stabilize below 20%. 
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FIG. 12. — The SFR of galaxies as a function of density for three different extinction corrections. The dashed lines are for the median and 75 '* percentile for 
the standard Kennicutt (1998b) formalism, assuming 1 magnitude of extinction. The solid lines are the median and 75'* percentile for the SFR-dependent dust 
correction using Hopkins et al. (2001). For the median, these two SFR indicators are very similar because the empirical corrections of Hopkins et al. (2001) are only 
significant for high SFR galaxies as shown by their difference in the 75'* percentile. The dashed-dot line shows the results of the Chariot & Longhetti (2001) model 
that uses information from six emission lines to constrain the dust properties of the galaxies. 



